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^ ■ Abstract 

■^ ' We discuss the stability of strangelets — quark droplets with strangeness — in 

the Nanibu-Jona-Lasinio model supplemented by a boundary condition for quark 
confinement. Effects of dynamical chiral symmetry breaking are considered properly 
inside quark droplets of arbitrary baryon number. We obtain the energy per baryon 

^ ■ number of quark droplets with baryon number from one to thousands. It is shown 

C<^ '. that strangelets are not the ground states as compared with nuclei, though they can 

^ I be locally stable. 

o . 

^ ■ 1 Introduction 

o 

'^ i The strange matter, which is a quark matter composed of the u, d and s quarks, has been 

I \ one of the interesting subjects in the quark and hadron physics [H El Ej- This subject 

Q^' is also related with other fields, such as the topics of compact stars, the early universe, 

T^ • heavy ion collision phenomena and so on. 

However, the physics of the strange matter has been a difficult problem of QCD. Ef- 
KA \ fective models of QCD were used from the early stage. The MIT bag model has been 

J^ \ often used to show the stabihty of the strange matter Pfl |21 El El- In this model, it is 

assumed that the strange matter is approximated as almost free quark gas with current 
quark masses. There, it was shown that the strange matter can be stable as compared 
with the ud quark matter. This is because the s quark current mass is small as compared 
with the Fermi energy of the ud quark. Moreover, it was argued that the strange mat- 
ter could be more stable than nuclear matter. Quark droplets with finite volume with 
strangeness, which we call strangelets, could be also more stable than atomic nuclei. If the 
strangelets are stable particles, they should be observed as almost neutral heavy particles 
in experiments. Indeed, several candidates of exotic particles have been observed in the 
cosmic rays 13 El El UHl HI]- There are also experiments to search for such exotic 
particles in heavy ion collision experiments fU ^1 E] . 

In the arguments in pi El El E] , it was assumed that chiral symmetry is restored in 
the strange matter. However, when we allow chiral symmetry breaking, the stability of 



cd 



the strange matter is not necessarily concluded. This is because s quarks acquire a large 
dynamical mass by chiral symmetry breaking. In P^ IT^ ITBj . the stability of the bulk 
strange matter was discussed by considering dynamical breaking of chiral symmetry in 
the Nambu-Jona-Lasinio (NJL) model, in which the quark dynamical mass depended on 
the baryon number density. In these discussions, however, it was shown that the strange 
matter could not be absolutely stable. The main reason of the instability of the strange 
matter is that the dynamical mass of the s quarks is larger than the ud quark Fermi 
energy. It was pointed out that the dynamics of chiral symmetry breaking is related 
with the formation of quark droplets [Hj . The stability of strangelets were discussed in a 
constituent quark model with one-gluon-exchange interaction [iSlCniEn]- There, it was 
assumed that the quarks have constant constituent mass. Thus, chiral symmetry breaking 
is an important ingredient to discuss the stability of the strange matter. 

It was also pointed out that chiral symmetry breaking plays an important role, not 
only in the bulk quark matter, but also in strangelets of finite volume \21\ |22l IIB] • In 
[221 I2n] , it was argued that the chiral symmetry breaking in the strangelets depended on 
the size of the strangelet by using the NJL model supplemented by the quark confinement. 
In this argument, it was assumed that the quarks interacting by the NJL interaction are 
confined in a MIT bag as a strangelet. This model is called the NJL+MIT bag model. 
In our previous work J22t l23] . instead of solving the boundary condition of the MIT 
bag exactly, the Multiple Reflection Expansion (MRE) [21] was used as an approximate 
method with an expectation that the strangelets have large baryon number. There, the 
surface effects were included by calculating the density of states of quarks modified by 
the boundary condition. It was shown that the chiral symmetry of the s quark sector in 
the strangelet is restored by the surface effect up to the baryon number of thousands, and 
strangelets could be more stable than the ud quark droplets. However, it was not possible 
to compare the stability of the strangelets with baryons (nucleons, hyperons and so on), 
since the mass of a single baryon was not calculated due to the limitation of the MRE 
approximation. 

In this paper, we investigate the NJL+ MIT bag model as one of the effective models 
of the baryon and the quark droplet. We consider directly the discrete quark energy 
levels by solving the boundary condition, which implements the quark confinement. The 
advantage of the present method is that we can construct a solution for all baryon number 
A including the nucleon and hyperon. This is contrasted with the previous discussion by 
the MRE which can be applied to the case of relatively large A J2S1- In the present study, 
we discuss the absolute stability of the strangelets as compared with baryons. In addition, 
we check the validity of the approximation of the MRE. 

Contents of this paper are the foUowings. In Section 2, we formulate the NJL + MIT 
bag model to describe chiral symmetry breaking in a strangelet. In Section 3, we discuss 
our numerical results by solving the discrete quark energy levels. There, the efficiency of 
the MRE approximation is investigated. We discuss the masses of a single baryon and 
quark droplet, and compare the stability between them. In Section 4, we conclude our 
discussion. 



2 Formulation 



The dynamical breaking of chiral symmetry is one of the important aspects of the low 
energy QCD. As an effective theory of QCD, we consider a four-quark interaction of the 
N JL model for chiral symmetry breaking j23] . For a quark droplet of finite size, valence 
quarks are assumed to exist inside the quark droplet. In order to implement this condition, 
we assume that quarks are confined in the MIT bag. Namely, we impose the MIT bag 
boundary condition at the bag surface of the quark droplet. With these assumptions, we 
use the lagrangian 
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e{R-r)-]-il)i)5{r-R), (1) 



where V = {u,d,sy is the quark field, and m° = diag(m°, mj], m°) current quark mass 
matrix. The second term in the first bracket is the NJL point-like interaction term 
invariant under U{Nf)L x [/{Nf)^, in which A*^ (a = 0, ■ ■ ■ , 8) are the Gell-Mann matrices 
normalized by trA'^A'' = 26"'''. We do not consider the 't Hooft term for explicit U{1)a 
breaking. In Eq. (P), the step function is multiplied in order to have quarks only inside 
the bag, and the last term realizes the MIT bag boundary condition for vanishing vector 
current flow at the bag surface |2El I2Z1 I2H1 ISHj- Formally, our model violates chiral 
symmetry at the bag surface. However, it is expected to be a good approximation that 
we do not consider the pion and kaon cloud for relatively large quark droplets. Details 
will be discussed in the later sections. Here, r is a distance from the center of the bag. It 
is assumed that the strangelet has a spherical shape with a bag radius R. We call Eq. (^ 
as the NJL+MIT bag model. In our model setting, the outside of the bag is described by 
the usual NJL model, which gives the chirally broken phase as vacuum. 

Chiral symmetry breaking is caused by the non-perturbative dynamics of the NJL 
interaction, in which qq correlation takes a non-zero expectation value. To investigate the 
qq condensate inside the bag, we use a mean field approximation for the scalar channel 
by replacing {qqY — > 2qq{qq) — {qq)"^ for q = u,d and s, respectively. Here (qq) is the 
expectation value of qq. Then, the lagrangian is approximated as 
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where we have defined 
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as an order parameter for chiral symmetry breaking, which also represents the constituent 
quark mass. We also have defined the dynamical mass of a quark by 
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--m +(j)g. 



(4) 



In this approximation, we obtain a single particle equation inside the bag (r < R) 

{i^-mg)q = 0, (5) 



with the boundary condition at the bag surface (r = R) 

— in ■ ^q = q, (6) 

which is obtained from the last term in Eq. (J2I). Here, 7 is spatial part of the gamma 
matrices. 

We solve Eq. (0) by using spherical basis set. The quark eigenstates of momentum 
p=\'p\, total angular momentum j = / + e/2 (e = ±l) and parity P=(— 1)^ are given by, 

^Uf^j^) =^ ( .^tU-n ) yu-)^ m 

where n is a unit vector with angular direction n = x/r and r = |x|. In Eq. (0), ji is 
the spherical Bessel function, and A/" is a normalization constant. The spinor harmonics 
yj^(n), which are eigenstates of the total angular momentum j = l+e/2 (e = ±l), is given 
by 

yjmin)=J2(^^^o^-^^\3^)'^U^)Xm-f„ (8) 

where V/^(^) are spherical harmonics and Xfi two component spinors. We write a quark 



energy E = Jp^+m? with quark mass niq. 

The boundary condition Eq. © of this quark wave function gives the momentum p in 
positive energy state by 

3i{pR)=e-^3i-H{vR). (9) 

and the momentum p in negative energy state by 

j^[pR) = -e-^ji{pR). (10) 

The absolute energy value of the state j^ in the vacuum is as same as that of the state j~^ 
in the valence. By solving the boundary condition (jU)) and ()10|). we obtain a sequence of 
states of a momenta p = PjPn{R, ""^q) with node quantum numbers ra as a function of the 
bag radius R and the quark mass rriq. The momentum pjP^iR, frtq) must be self-consistent 
with the mean field approximation Eqs. (jH)) and P|). 

We obtain the energy density e of strangelets with radius R by taking a sum of the 
energy over all states including positive and negative energy states in the bag. 
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where u = N^pin xA^^ = 2x3 = 6 is the number of degeneracy, and for positive energy state 



E,p,n{<P,) = ^|pf~HK + ^^^ (12) 
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and for negative energy state j^, 



EjP,n{.(t)q) 



P'^iP.n + ("^0 + 



(13) 



Here, PjP^n ^'Hd PjP^n ^I's the momenta in the valence and vacuum parts, which are given by 
the boundary condition @ and pO|). respectively. The Fermi momenta piT'g are obtained 
by 



PFq 



E=N. 
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for a given quark number Nq for each flavor q = u,d and s. Since the NJL model is a 
cutoff theory, we have introduced a regularization function g{p/A) for a sum over negative 
energy states. The results may depend on the choice of g{p/A). In the present work, we 
use a smooth function of the Lorentzian type 



^(P/A) 
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with a diffuseness parameter a. We consider that the outside of the bag is chirally broken 
phase descried by the NJL model. In Eq. (fTT|) . the energy density of the strangelet e is 
measured from the vacuum of the chirally broken phase in the absence of the bag. In the 
NJL model, the latter energy density is calculated by 
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is the value of the order parameter in the chirally broken 



where 0o = (0«o, 0do, 0so) 
vacuum. 

For later discussions, we write the energy density of the strangelet as a sum of the 
kinetic energy part of the valence quarks and the rest. 



5(0) 



VFq 



Y: E^F,n(0.)+5(0)- 
q=u,d,s jP^n 



(17) 



Here, the effective bag constant B{(j)) is the energy density of the vacuum inside the bag 
measured from Eq, 

B{<P) = eZM)-^o, (18) 

where we define the vacuum energy inside the bag 
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We mention that the effective bag constant B is not a constant independent of the size 
of the bag, but is a quantity depending on the mass of the quarks nig, or 0^, and the bag 
radius R. 

By using the energy density (fTT|) . the total energy of a strangelet is given by 

E = €V-a/R, (20) 

with the volume V = ^R^. The second term of this equation expresses the correction 
from zero point energy, center-of-mass motion and other effects ^U\ |M] . Then, we obtain 
a locally stable state of a quark droplet for a given baryon number A and strangeness 
S {S = —Ng and Nu = Nd). There, the order parameter 0^ and the bag radius R are 
determined by stationary condition of the energy E with respect to (pq (0„ = (j)d) and R, 
respectively; 

dE dE 

We have determined the model parameters as follows; the three dimensional momen- 
tum cut-off A = 0.60 GeV, the coupling constant GA^ = 4.61, the current quark masses 
m^ = r7i° = 0.0059 GeV and m° = 0.134 GeV. This parameter set reproduces the pion 
decay constant /7r = 0.093 GeV, pion mass m.,, = 0.139 GeV, kaon mass mx = 0.493 GeV, 
and dynamical quark mass rriu = rrid = 0.378 GeV in the chirally broken vacuum. The 
diffuseness parameter a in the regularization ()15|) is treated as a free parameter by fitting 
the baryon mass. We use the value a = 22 in Subsections 3.2 and 3.3. In Subsection 3.3, 
we discuss the a dependence of the baryon mass. 

3 Numerical results 

In this section, numerical results are investigated for the strangelets with the baryon 
number A. 

3.1 Comparison with MRE 

First, we compare the present result with the previous ones obtained by the Multiple 
Reflection Expansion (MRE). In the MRE, the surface term in the quark droplet is intro- 
duced as a correction to the bulk quark matter |21]. The MRE was often used for analysis 
of quark droplets in the MIT bag model |lj . 

In our previous work in [21], the MRE was applied to the strangelets with baryon 
number A > 10. Now, we can check the efficiency of the MRE approximation, by com- 
paring the number of valence quarks in the MRE with that obtained by calculating the 
discrete levels. In the bag model, we solve the boundary condition and (fTIH) . Then, 
for a bag with a radius R, the total fermion number in the momentum space up to the 
momentum A is obtained as 

N,,g{A) = Y.9{A-pjP^^). (22) 



On the other hand, the MRE gives 

NuREiA) = 2 1^ ^PMREip, m, R), (23) 

where the factor 2 is the number of degeneracy of spin and Pmre{p, i^, R) is the density 
of state of the MRE jH IM|; which is given as 

evr^ 127r2 

Pmre{p, m,R) = 1 + —fsip/m) + -r^fcip/^) + ■■■, (24) 



with 



1 / 2 
fs{x) = -^ 1 arctanx), (25) 
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We show the number of states Nhag{A) and Nmre{-A) for the fermion mass m = 0.3 
GeV for the bag radii R = 1 and 2 fm, in Fig. ^a) and (b), respectively. For comparison, 
we show also the number of states in the bulk matter 

iV,„,,(A) = 2 r^. (26) 
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The behavior like step function is due to the discreteness of the quark levels. In both fig- 
ures (a) and (b), we see that the result of the MRE agrees well with the averaged behavior 
of the present discrete levels. Furthermore, it is clear that the MRE approximation agrees 
better, as the bag radius increases. Thus, we can use the MRE approximation instead 
of the discrete level calculation for larger bag radii. In later discussions, we will confirm 
that the energy per baryon number of strangelets calculated from the discrete levels is 
consistent with that in the MRE approximation for large bag radii i? > 2 f m (see Figs. |21 
and|3|). 

3.2 Vacuum structure in the bag 

The boundary conditions Q and (fTUj) modify the quark energy levels in the bag, which 
is known to be an origin of the Casimir effect. This affects not only the valence quarks 
with positive energy, but also the vacuum quarks with negative energy. Accordingly, the 
properties of the dynamical chiral symmetry breaking is also changed. In this subsection, 
we discuss the vacuum structure and the dynamics of chiral symmetry breaking in the 
bag . 

To start with, we consider an empty bag with only vacuum quarks. The energy of the 
vacuum inside the bag is given by the effective bag constant -8(0) of (fT8|l . This quantity 
is obtained by taking a variation with respect to the order parameter cj) = {ipu, (f>d, 4>s) 
for the bag with a given bag radius R. The resulting effective bag constant B and the 
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Figure 1: The number of states of quarks {Nhag (solid lines), Nmre (dashed lines) and N^i^ik 
(dashed-dotted lines)) in the bag as functions of the three dimensional momentum cut-off A for 



radii (a) R 
GeV. 



1 fm and (b) i? = 2 fm , respectively. The fermion mass is set to be ruq = 0.3 



order parameter of chiral symmetry breaking in s quark sector </)s are shown in Fig. |21 (a) 
and (b), respectively, as functions of the bag radius R. The order parameter in the ud 
sector is zero for the region of the bag radii -R < 6 fm. We compare the results obtained 
by calculating the discrete levels (solid lines) and by the MRE (dashed lines). For later 
convenience, we also show the effective bag constants and the radii of quark droplets filled 
with valence quarks by blobs. The corresponding baryon numbers are shown in the same 
figure. These values are obtained by taking a variation of the energy with respect to the 
quark mass and the bag radius. The quark droplets with valence quarks are discussed in 
Section 3.3. 

In Fig.Efa), it is shown that the effective bag constant is a function of the bag radius R. 
As a bulk property, the effective bag constant becomes small as the bag radius increases. 
Indeed, we see a good agreement between the MRE and the present calculation with the 
discrete levels for larger bag radii R>2 fm. The state in the bag in the limit of the large 
bag radius i? ^ oo is chirally broken bulk vacuum, because the surface effect of the bag 
can be neglected. In this limit, we obtain B = form Eq. ^T^ by substituting = 0o- 

On the other hand, in the limit of small bag radius i? — *> 0, the effective bag constant 
becomes a positive quantity. Let us investigate slightly more of the vacuum structure in 
the small bag radius. The highest state in the negative energy sea in the bag is the l/2~ 
state with energy -E1/2- = — 2.04/_R. As the cut-off is introduced in the bag vacuum part 
in Eq. ()19j) . no level exists for bag radii R < 2.04/A(= 0.67 fm), since all the states are 



below the cut-off A. Therefore, we find slagi^') = J2 



(j>=u,d,: 



74(7 from Eq. (HH), which 



follows 0q = as chirally restored phase. Consequently, we obtain 

B = -6o (i?< 0.67fm), (27) 

which is a constant value independent of the bag radius. We show that —Eq is a positive 
quantity. For this purpose, we rewrite the vacuum energy Eq as 

Eo = -B,^ik + eZiki<P = 0), (28) 

where we define 

B,^ik = EZU<P = 0)-eo. (29) 

This is the energy density difference between the chirally broken and restored phases in 
the bulk space without considering the boundary condition. We also obtain the energy 
density of the vacuum quark in the chirally restored phase. 



^vac 



(2^ 



eiZIki^ = 0) = -u ^^p2 + m%g(p/A). (30) 



It is clear that EZfi.((j) = 0) is a negative quantity, and Bt,uik is also a positive quantity. 
Therefore, Eq is a negative quantity. Finally, we conclude that the effective bag constant 
B = —Eq is a positive quantity. 

The above analysis seems to suggest that the effective bag constant i? is a monotoni- 
cally decreasing function of the bag radius R. However, there is a sharp decrease of the 
effective bag constant around the bag radius i? ~ 0.7 fm in the discrete level calculation, 
which is not shown in the MRE, as shown in Fig. |2ta). This behavior is understood by 
considering the l/2~ state in the vacuum part in the bag. As shown in Eq. (jTTj) . we have 
a constant effective bag constant B = —Eq for the bag radii R < 2.04/A, because no 
quark state exists anymore in the negative energy of |pjp,„| > A. On the other hand, for 
the bag radius R > 2.04/A, several states, such as l/2~, 3/2+ and so on, are taken into 
account in the calculation of the effective bag constant. Especially, for the bag radius 
2.04/A < R< 3.20/A (0.67fm < R< l.OSfm), only the 1/2" state appears in Eq. (P|) . 
The next state 2/3"*" with energy £"3/2+ = — 3.20/i? does not exist in the vacuum in this 
bag. Therefore, for these bag radii, we obtain the effective bag constant 

B ~ - '^'^'^^'/^~ g{2M/AR) - e^ (0.67 f m < i? < 1.05 fm), (31) 

where A^i/2- = 18 is the number of degeneracy of l/2~ state including the m, d and s 
quarks, V = ^R^ the volume of the bag. Here, chiral symmetry of all flavors is restored 
as shown in Fig. |21^b), and quark masses are current ones. The energy of m, d and s quarks 
are almost degenerate, since the masses are smaller than the quark momenta inside the 
bag, p ~ 2/R ~ 0.5 GeV >> ml^^. Now, the effective bag constant (j^ is an increasing 
function B ~ —1/R'^ of the bag radius R. Consequently, from Eqs. (J27jl and (J3T| . we 
obtain the sharp decrease of the effective bag constant B around i? ~ 0.7 fm. 
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Figure 2: (a) The effective bag constant i? as a function of a radius of the vacuum bag. The blobs 
correspond to points of baryon number A as denoted by the labels, (b) The order parameter of 
chiral symmetry breaking in s quark sector (pg as a function of the bag radius R. In both figures, 
the solid lines are the present results of discrete levels, and the dashed lines are the results by 
the MRE method. 

3.3 Stability of strangelets 

In this subsection, we discuss the stability of strangelets with baryon number A and 
strangeness S. First, we show a global behavior of the quark droplets including the 
baryon. Next, we discuss the structure of quark droplets in detail. 



3.3.1 General properties of quark droplets 

First, we outline our result of the baryon number dependence of the properties of the 
quark droplets. In Fig. IHl we show the resulting energy per baryon number E/A as a 
function of baryon number A, which is obtained by taking a variation with respect to the 
chiral condensate (pu, 4>d and (ps, the bag radius R. We also consider the variation with 
respect to the strangeness S = —Ng. The solid line is obtained by considering the discrete 
energy levels, and the dashed line is the previous result obtained by the MRE method 
[2H] • Here, we define the strangeness fraction r^ as a ratio of the s quark number to the 
total quark number, Vg = Ng/SA. 

From Fig. IHl in the large baryon number A > 10, the present result of the discrete 
levels agrees well with that in the MRE. This justifies the validity of using the MRE 
method for the baryon number A > 10 J2S1- On the other hand, in the small baryon 
number A < 10, the MRE does not agree with the present calculation. The discrete level 
calculation reproduces nucleon mass rriN = 1.10 GeV at A = 1. Instead, the mass of 
the quark droplet in the MRE becomes an unphysical value above 2 GeV as the baryon 
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Figure 3: The energy per baryon number E/A as a function of baryon number A. Tlie solid 
line is by the discrete energy level, and the dashed line is by the MRE method. 

number approaches A = 1. 

We show the resuhing strangeness fraction r^ as a function of the baryon number A in 
Fig. m Here, we define the strangeness fraction r^ as a ratio of the s quark number to the 
total quark number, r^ = Ns/3A. It is shown that the quark droplets with baryon number 
5 < A < 10'^ have the finite strangeness fraction r^ ~ 0.3. Therefore, the strangelets can 
exist as local minima for the limited baryon number 5 < A < 10^. On the other hand, the 
ground state of the quark droplets with baryon number A > 10^ is the ud quark matter. 
The value of the critical baryon number A ~ 10^ was also obtained in the calculation by 
the MRE 1^. The quark droplets in the limit of the large baryon number is the bulk 
quark matter of ud quarks. This is consistent with the results obtained by the NJL model 

[miiniiiE]. 

The chiral symmetry breaking of the s quarks affects the stability of the strangelets. 
The dynamical quark mass rriu and m^ of the ud and s quark are obtained as functions 
of the baryon number A in Fig. El The ud quark Fermi energy epu is shown in the same 
figure. The chiral symmetry of the ud and s quarks is restored for A < 10^. The s quarks 
are contained in the ground state of the quark droplets, since the relation 



rris < epu, 



(32) 



is satisfied. On the other hand, the chiral symmetry of the s quarks is broken at A ~ 10^, 
where the dynamical s quark mass is equal to the ud quark Fermi energy. Therefore, s 
quarks do not appear in the quark droplets with baryon number A > 10^. 

Now, let us discuss the structure of the quark droplets for three ranges of the baryon 
number. 
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Figure 5: The ud quark Fermi energy epu and the s quark dynamical mass mg as functions of 
the baryon number A. 
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Figure 6: (a) The energies E oi a, baryon {A = 1) for the strangeness S" = 0, —1, —2 and —3 as 
functions of the bag radius R. (b) The energy levels of eigenstate J^ of the u, d (left row) and s 
(right row) quarks in the valence (upper) and the vacuum (lower) part in the bag of the baryon 
with strangeness 5 = and radius R = 0.75 fm. The numbers in the circles are the number of 
the quarks occupying the levels. 

3.3.2 The case of 1 < A < 4 

We first consider a quark droplet with A = 1, namely a baryon. In Fig. IHl^a), we show the 
energies E for baryons with each strangeness S' = 0, —1, —2 and —3 as functions of the 
bag radius R. There are local minima E = 1.10, 1.16, 1.23 and 1.30 GeV at bag radius 
R = 0.75 fm for S = 0, —1, —2 and —3, respectively. The ground state is the nucleon with 
the mass rriN = 1-10 GeV and the bag radius R = 0.75 fm. The inside of the baryon is 
the chirally restored phase (see Fig. |3J. 

Now, we discuss the quark level structure in the baryon shown in Fig. IHfb), where we 
show the energies of the quark eigenstates in units of GeV on the left of the levels and 
the corresponding total angular momentum j and parity P on the right. The numbers 
in the circles on the levels are the number of quarks that occupy the levels. The dashed 
lines in the vacuum part indicate the momentum cut-off in the regularization. The states 
below the cut-off are hidden in the shadow. 

The three valence quarks occupy the state 1/2+ in the positive energy part of the ud 
quark sector. The state 1/2+ of the s quark sector is not occupied, since its energy is 
larger than that of the ud quark due to the difference of the current masses. Therefore, 
the ground state of the baryon is non-strangeness state of 5* = 0. Concerning the vacuum 
quarks, only the l/2~ state is included within the momentum cut-off as discussed in 
Eq. ()3H1 . The next level of 3/2"*" and more states existing below are not inside the cut-off 
A. 

We stress that the shell structure in the vacuum part inside the bag is important for 
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hadrons 


S 


E (cal.) [GeV] 


E (exp.) [GeV] 


A^ 





1.10 


0.939 


A 





1.10 


1.230 


A 


-1 


1.16 


1.115 


S 


-1 


1.16 


1.190 


z: 


-2 


1.23 


1.314 


Q 


-3 


1.30 


1.672 



Table 1: The energy of baryon with strangeness S". We denote i?(cal.) and£'(exp.) as calculated 
and experimental values, respectively. 

the baryon mass. Indeed, the baryon mass was not obtained by the MRE method in |2S!, 
because the MRE method was not appropriate to discuss the properties of a small size 
quark droplet due to the limitation of the approximation of this method. In the present 
method, despite the simphcity of the model, we have obtained the masses oi A = 1 
baryons, especially of the nucleon reasonably well. 

We have calculated masses of other SU(3) baryons. The result is shown in Table Q 
as compared with experimental values. We obtain E\ = E-^ = 1.16 GeV, E^ = 1.23 
GeV and Eq = 1.30 GeV. These masses are close to the experimental values. In our 
calculation, the octet and decuplet baryons, for instance N and A, are degenerate; spin 
dependent residual interaction is responsible for the octet-decuplet splitting. We check 
that the Gell-Mann-Okubo mass formula 



2(A^ + S) = 3A + S 



(33) 



holds within 1.1 % in our numerical result. 

Now, let us turn to the discussion for A > 2. Since the 1/2+ state can accommodate 
twelve quarks at most, the ground state of the quark droplets with the baryon number 
A = 2, 3 and 4 is still the non-strange state; the ud quark droplets. Then, we obtain 
the energy of the ud quark droplets given by local minima as functions of the bag radius, 
as we discussed in the case A = 1. Let us see the resulting energy per baryon number 
E/A of the quark droplets with baryon number A in Fig. El We obtain E/A = 1.36, 1.44 
and 1.48 GeV for A = 2, 3 and 4, respectively, which are larger than the nucleon mass 
E = lAO GeV. Therefore, the ud quark droplets of 2 < A < 4 are not stable as compared 
with the y4-nucleon system. The reason for the instability of the ud quark droplets is the 
increase of the effective bag constant B for the bag radius R for 0.7 f m < i? < 1.05 fm as 
shown in Fig. |2fa). For example, the quark droplets with baryon number A = 2, 3 and 4 
have the bag constant S^/^ = 0.077, 0.10 and 0.12 GeV, respectively, while the baryon of 
A = 1 has B^^^ = 0.036 GeV. Therefore, the energy per baryon number of the ud quark 
droplets becomes larger than the baryon mass. 
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3.3.3 The case of 5 < A < 10^ 

Next, let us discuss the quark droplets with baryon number A > 5. For these baryon 
numbers, we obtain the strangelets as the ground state of the quark droplets. Let us 
investigate the case oi A = 5. In Fig. E^a), we plot the energies per baryon number E/A 
of the quark droplets oi A = 5 for strangeness 5* = 0, —3 and —6 as functions of the 
radius R. At the local minima, we obtain E/A = 1.68, 1.54 and 1.58 GeV for strangeness 
S = 0,-3 and —6, respectively. The most stable state is the strangelet with the finite 
strangeness S = —3. The bag radius is R = 0.76 fm. We can understand this result 
of the finite strangeness by considering the structure of the valence quark levels in the 
bag as shown in Fig. [7|^b). Here, we use the current quark masses to obtain the energy 
levels, since chiral symmetry is still restored at A = 5 (see Fig. EI). The twelve quarks 
occupy fully the lowest state 1/2+ of the ud quark sector. We may put the remaining 
three quarks in the next level 3/2~ in ud quark sector. However, the energy -E'3/2- = 0.83 
GeV of the ud quark sector is larger than the energy -E1/2+ = 0.60 GeV in the s quark 
sector. Therefore, the state 1/2+ of the s quark is occupied by the remaining three quarks. 
Therefore, the ground state of the quark droplet with baryon number A = 5 is strangelets 
with strangeness S = —3. 

We have investigated the energy of the strangelet for larger A, and found that the 
strangelets can be the ground state as quark droplets with A < 10^. From Fig. 13 the 
energies per baryon number E/A of the strangelets of 5 < A < 10^ are 1.5 GeV < E/A < 
1.8 GeV, which are larger than the nucleons mass E^^f = 1.10 GeV and the Lambda mass 
E/^ = 1.16 GeV. Therefore, the strangelets are not stable against normal nuclei and 
hypernuclei. 

3.3.4 The case of A > 10^ 

For the baryon number A > 10^, the ground state of the quark droplet is not strangelet, 
but the ud quark droplet. This is because the chiral symmetry in the s quark sector is 
broken. As shown in Fig. the dynamical mass rris of the s quark becomes larger than 
the ud quark Fermi energy epu for A > 10^. There, the condition for the weak decay (jH^ 
is not fulfilled. Therefore, the s quarks cannot appear in the ground state of the quark 
droplets. 

From the discussion above, it is shown that the baryon number of stable strangelets 
is limited in 5 < A < 10^. This conclusion is qualitatively different from the result in the 
MIT bag model; all strangelets of the baryon number, say A ^ 10, are the most stable 
state as compared with nuclei jSlll]. In contrast, our result by the NJL+MIT bag model 
indicates that the strangelet can exist as a quasi-stable state rather than the true ground 
state. 

3.3.5 Choice of the parameter a 

Up to now, we have used a parameter set reproducing the properties of the infinite volume 
vacuum. However, the values of the diffuseness parameter a in the Lorentzian regulariza- 
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Figure 7: (a) The energies per baryon number E/A of quark droplets with baryon number 
A = 5 for the strangeness S* = 0, —3 and —6 as functions of the bag radius R. (b) The energy 
level of eigenstate J^ of the n, d (left row) and s (right row) quarks in the valence (upper) and 
the vacuum (lower) part in the bag of baryon number ^ = 5, strangeness S = —3 and radius 
R = 0.76 fm. The number in the circle is the number of the occupying quarks in each state. 

tion (fT3j) cannot be determined only from the properties of the infinite volume system, 
because such a system is insensitive to the parameter a. In order to determine a, let us 
look at the a-dependence of the nucleon mass and radius. 

We show nucleon mass E, the bag radius R and the effective bag constant B for several 
diffuseness parameter a in Table El As the diffuseness parameter changes from a = 20 
(smooth) to 25 (sharp), the baryon mass E and the effective bag constant B^/^ changes 
from E = 1.32 GeV to 0.83 GeV and from B^/^ = 0.177 GeV to -0.186 GeV. For a = 22, 
we find ii^ = 1.10 GeV which is the average value of the masses of the nucleon and the 
delta. This is the reason why we have used a = 22 in the preceding discussions. On the 
other hand, the nucleon radius remains around i? ~ 0.75 fm for the range of a shown in 
Table El We also checked that the Gaussian regularization with an appropriate parameter 
set. However, the nucleon mass was not reproduced in the Gaussian type. 

The energy of the quark droplet in the present framework depends on the parameter a 
particularly for small A region. However, the present result for A > 10 is rather parameter 
independent, and the result oi A = 1 should be like the one of the nucleons and the delta, 
as it was reproduced in the present treatment by suitably choosing a. Therefore, we 
expect a smooth interpolation of the two points, for instance {A = 10, E/A = 1.8 GeV) 
and {A = 1, E = 1.10 GeV) of Fig. |21 although we need further refined treatment for the 
detailed behavior for small A region. 
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a 


nucleoli 


ud quark matter 


E GeV] 


R [fm 


5^/4 GeV 


e/ub GeV 


ns/nBo 


5^4 [GeV] 


20 


1.32 


0.747 


0.177 


1.121 


2.6 


0.174 


21 


1.20 


0.748 


0.147 


1.117 


2.6 


0.174 


22 


1.10 


0.749 


0.0362 


1.114 


2.6 


0.173 


23 


1.00 


0.749 


-0.144 


1.112 


2.5 


0.173 


24 


0.90 


0.749 


-0.167 


1.110 


2.5 


0.173 


25 


0.83 


0.750 


-0.186 


1.107 


2.5 


0.173 



Table 2: We show nucleoli mass E, the radius R and the effective bag constant B oi a. single 
baryon for several diffuseness parameter a in the left side. In the right side, the energy per 
baryon, the stable density and the effective bag constant of ud quark matter are shown. 

3.3.6 Other quantities 

We discuss radii and baryon densities of strangelets. In Fig. we show the bag radius 
R of the strangelets as a function of the baryon number A. The solid line shows the bag 
radius which is obtained by the discrete level calculation. We fit our numerical result 
by an approximate relation R = roA^^^ with a constant tq, in the dashed line. Here, we 
obtain the constant vq = 0.56 fm by fitting our result. Then, the averaged baryon number 
density is ub = A/V ~ Subo, where ubq = 0.17 fm~^ is the normal nuclear matter density. 
Hence, the strangelets are more compact objects as compared with ordinary nuclei. 

A strangelet existing as a locally stable state has a small electric charge as compared 
with the normal nuclei. The electric charge is given hj Q = {2/3)Nu~{l/3)N(i—{l/3)Ns = 
0.5(1 — 3rs)A. From Fig. 01 we obtain the charges per baryon number < Q/A < 0.3 for 
the strangelets of 5 < A < 10^, which are smaller than the values in the normal nuclei 
Q/A ~ 0.5. Therefore, the strangelets would be observed as heavy particles with small 
electric charges. 

3.4 Bulk quark matter 

In this subsection, we discuss the stability of the bulk quark matter in comparison with 
the nuclear matter. In |14[ IT^ IT^ , it was shown by using the N JL model that the ground 
state of the quark matter is not the strange matter, but the ud quark matter. This result 
agrees with our present result that the quark droplets with large baryon number A > 10'^ 
are the ud quark droplets, not the strangelets. In [U UHl HE] , however, it was also not 
possible to compare the stability of the ud quark matter with that in the nuclear matter, 
since the baryon mass was not derived in the NJL model in their treatment. In contrast, 
our present model can reproduce the baryon mass, which can be compared with the energy 
of the bulk quark matter. 

We consider that the bulk quark matter in the present model (0). This is obtained 
either by neglecting confinement term in Eq. (^, or by taking the limit of the infinite 
radius. Then, in the mean field approximation Q and Q, the energy density of the bulk 
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Figure 8: The solid line is the radius K of strangelets as a function of baryon number A. The 
dashed line represents R = vqA^'^ with tq = 0.56 fm. 



quark matter is given by [HI [T31 Hn] 



q=u,d,! 






{27ry 



The quark number density is 



n„ 



9 (P/ A) ^P'^ + (mog + 0g)^ 






£0- 



(34) 



(35) 



for each flavor q = u, d and s. In Eq. (I34|l . instead of the sharp cut-off, we have a smooth 
cut-off by g{p/A). 

In Fig. Efa), we show the energy per baryon number e/ub as functions of the baryon 
number density ub = {nu + na + ns)/3 for the diffuseness parameters a = 20, 22 and 24. 
We also show the strangeness fraction r^ obtained by a variation of the energy per baryon 
s/ub in Fig. IHl^b). The ground state of the quark matter is ud quark matter for lower 
density ub ^ lO.Snso, because the chiral symmetry of s quark sector is broken at the 
vacuum and the lower density. On the other hand, the ground state of the quark matter 
is strange matter for higher density ub ^ lO.buBo, since the chiral symmetry of s quark 
sector is restored. In fact, the strangeness fraction takes a finite value at the density 
Ub > lO.Snso- 

From Fig. El^a), there are stable states of the quark matter around around ub/ubo = 
2.5. This is the ud quark matter, not the strange matter. There, we obtain the energy 
per baryon number and the baryon number density of the stable ud quark matter. We 
summarize our results in Table |21 The energy per baryon is e/ub = 1.10 — 1.12 GeV, 
the stable density is ub/ubq = 2.5 — 2.6 , and the effective bag constant is B^^^ = 
0.173 — 0.174 GeV. These quantities in the bulk quark matter are not affected by the 
diffuseness parameter a so much as in the baryon. 
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Figure 9: (a) The energy per baryon number e/ub as functions of the baryon number density ub- 
(b) The strangeness fraction r^ as functions of the baryon number density hb- ubo = 0.17 fm~ 
is normal nuclear matter density, a is diffuseness parameter in the momentum regularization. 

Let us compare the energy E of baryon and the energy per baryon number e/ub in the 
ud quark matter. For the diffuseness parameter a > 22, the ud quark matter is unstable 
against the baryon. On the other hand, for the diffuseness parameter a < 21, the baryon 
mass E exceeds the energy per baryon number e/hb of the ud quark matter. Therefore, 
for the diffuseness parameter a < 21, we can expect that the quark matter can be stable 
state rather than the nuclear matter. From Fig. El the ud quark droplets in the limit of 
the large baryon number can be more stable than the baryon system, when the energy 
per baryon of the ud quark matter is smaller than the baryon mass. Even in that case, 
however, the strangelets are not absolutely stable objects due to the large energy per 
baryon as compared with the baryon mass. 

4 Conclusion 

In this paper, we discussed the stability of the quark droplet by considering dynamical 
chiral symmetry breaking in a finite volume region of the droplet. We considered chiral 
symmetry breaking by using the NJL model supplemented by the MIT bag boundary 
condition for quark confinement. The boundary condition prepared the basis set in the 
bag. Then, we obtained the energy per baryon number of the quark droplets by solving 
self-consistently the equation of state in the mean field approximation. Masses and radii 
of various baryons were also discussed in this model. 

We obtained the result that the strangelets are stable as quark droplets for the baryon 
number 5 < A < 10^. This is because the chiral symmetry of the s quark sector is restored 
for these baryon numbers. We also obtained the baryon mass Ej^f = 1.10 GeV and the 
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radius i? = 0.75 fm as a quark droplet with baryon number A = 1, which was identified 
with the nucleon. However, the energy per baryon number of strangelets are larger than 
the baryon mass. Therefore, we conclude that the strangelets can not be the ground 
state in the present framework. We also discussed the stability of the bulk quark matter 
against the nuclear matter. The baryon mass depends on the diffuseness parameter in the 
Lorentzian regulator, while the energy of the bulk matter did not. By comparing baryon 
mass and the energy per baryon number in the stable ud quark matter, it was concluded 
that the ud quark matter could be more stable than the baryon. 

We also confirmed that the MRE method reproduces the results of the present cal- 
culation using discrete levels directly for the baryon number A > 10. This justifies the 
use of the MRE for larger A, while for smaller A the shell structure of the discrete levels 
becomes important. 

We have still some possibilities that the energy of strangelets is modified by other 
effects which were not fully considered in our discussion. As one of the most important 
terms, we need to consider the 't Hooft term for explicit U{1)a breaking. We may need to 
consider the meson clouds around the bag also. In fact, from the view of chiral symmetry, 
the TT and K mesons should be introduced at the bag surface in order to recover chiral 
symmetry broken at the bag surface. In the present work, we neglected them. This 
approximation is expected to be valid for systems with large baryon number. For instance, 
in the previous studies of the chiral bag model, it was shown that the pion effect for 
A = 2 is significantly smaller than the nucleon {A = 1) |321 inBj- Therefore, at least 
for A > 2, we expect that the present result will not be affected so much by the meson 
cloud. Nevertheless, the inclusion of the meson cloud in the model is very important 
in particular for the interpolation of the present analysis to the small A region (1 < 
A < 10). We are currently working out such a problem. Another interesting effect is the 
diquark correlation in the quark droplets [SH ESI IHE] • In recent studies of high density 
matter, it has been discussed that the true ground state of the quark matter may be color 
superconductivity state. Such effects would be necessary to further development of the 
study of the strangelets. 
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